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1 Introduction 

A number of possible mechanisms have been suggested to generate density inhomo¬ 
geneities in the early Universe which could survive until the onset of primordial nu¬ 
cleosynthesis and affect the abundances of the isotopes produced in the big bang [1]. 
The possibility that the changes in the abundance pattern might be drastic enough to 
reconcile a closure density of baryons with primordial nucleosynthesis constraints, or to 
produce a characteristic signature of a phase transition in the early universe has proven 
exciting enough to inspire a considerable amount of work on inhomogeneous nucleosyn¬ 
thesis. In this work we are not concerned with how the inhomogeneities were generated 
but we want to focus on the effect of such inhomogeneities on primordial nucleosynthesis. 
One of the proposed signatures of inhomogeneity, the synthesis of very heavy elements 
by neutron capture, is analyzed for varying baryon to photon ratios rj and length scales 
L. A detailed discussion is published in [2]. Preliminary results can be found in [3]. 


2 Method 

After weak decoupling the vastly different mean free paths of protons and neutrons 
create a very proton rich environment in the initially high density regions, whereas the 
low density regions are almost entirely hlled with diffused neutrons. Since the aim of the 
present investigation was to explore the production of heavy elements, we considered 
only the neutron rich low density zones. High density, proton rich, environments might 
produce some intermediate elements via the triple-alpha-reaction, but will in no case 
be able to produce heavy elements beyond iron. However, we included the effects of 
the (back) diffusion of neutrons into the proton rich zones. Using a similar approach as 
introduced in [4, 5], the neutron diffusive loss rate /t is given by 

^ + 0.716T9)V25-i (1) 

{a/ajcmMeV 

in the temperature range 0.2 < Tg < 1. Thus, the only open parameter in the neutron 
loss due to diffusion is the comoving length scale of inhomogeneities {d/a). Small sep¬ 
aration lengths between high density zones make the neutron leakage out of the small 



low density zones most effective. Large separation lengths make the neutron leakage 
negligible. (For a detailed derivation of Eq.(l), see also [2]). 

Our reaction network consisted of two parts, one part for light and intermediate nuclei, 
the second part being an r-process code. For light and intermediate nuclei from neutrons 
and protons to krypton (Z=36), from stability to the neutron drip line, we made use of 
a general nuclear network (of 655 nuclei) which includes neutron, charged particle, and 
photon induced reaction as well as weak reactions. (For details of the included rates see 
Appendix B and Tables 1, 2, Al, and A2 in [2]). 

The second part was an r-process code that determines the abundances of heavy nuclei. 
This network extends up to Z = 114 and contains all (6033) nuclei from the so-called 
valley of beta-stability to the neutron-drip line (see also [6]). The neutron capture rates 
were calculated with statistical model methods and the beta decay rates were taken 
from [8], where experimental values were not available (see also [7]). For the calculations 
in this paper we also introduced (beta delayed) hssion of heavy nuclei, as calculated by 
Thielemann, Metzinger, and Klapdor (1983) [9] with the Howard and Moller (1980) [10] 
hssion barriers and masses (see [7]). These two networks were coupled together such that 
they both ran simultaneously at each time step, and the number of neutrons produced 
and captured was transmitted back and forth between them. 


3 Results and Discussion 

The choice of an initial neutron abundance of = 1 (i.e. only neutrons, which is the 
most favorable condition for the formation of heavy elements) in the low density region 
leads to a density ratio piow/pb = 1/8 [2]. This leaves as open parameters the baryon to 
photon ratio rj = rib/n^ = 10“^°pio and the comoving length scale {d/a). Four sets of 
calculations have been performed, employing pio values of 416, 104, 52, and 10.4. Using 
the relation [2] 

ilbhlo = 1.54 X 10-2(T,o/2.78iF)Vo , (2) 

with the well known present temperature of the microwave background and the 
Hubble constant Ho = h^o ^ 50 km s“^ Mpc“^, this corresponds to possible choices of 
{h^Q,flb) being (2.5,1), (1.3,1), (1,0.8), and (1,0.16). The range covered in pio extends 
from roughly a factor of 2.2 below the lower limit to a factor of 13 above the upper limit 
for f] in the standard big bang. For each of the p-values we considered four different 
cases of d/a: (0) oo, resulting in negligible neutron back diffusion, (1) lO^'® cm MeV, (2) 
10®'^ cm MeV, and (3) 10^'^ cm MeV. (This corresponds to distances between nucleation 
sites of oo, 2700, 270, and 27 m, respectively, at the time of the quark-hadron phase 
transition). The resulting abundances for heavy elements are shown in Tab. 1. One 
notices the exponential increase in r-process abundances with increasing rj. This is due 
to “hssion cycling”, whereby each of the hssion fragments can capture neutrons and 
hnally form again heavy nuclei, which are also prone to hssion [11]. This is of particular 
importance in environments with a long duration of high neutron densities, and was 
therefore suggested as relevant to primordial nucleosynthesis in neutron rich zones of an 
inhomogeneous big bang [4]. In contrast to the operation of the r-process in explosive 



stellar environments, confined to a few seconds, this process in the neutron rich regions 
associated with an inhomogeneous big bang is only limited by the neutron half-life and 
can go on for an extended period of time. One of the remarkable features of an r-process 
with hssion cycling is that the production of heavy nuclei is not limited to the r-process 
flow (neutron captures and beta-decays) coming from light nuclei, but requires only a 
small amount of hssionable nuclei to be produced initially. The total mass fraction of 
heavy nuclei is doubled with each hssion cycle and can thus be written as X^. = T^Xgeed- 
Here, n is the number of hssion cycles and Xgeed denotes the initial mass fraction of 
heavy nuclei. The ehectiveness of hssion cycling can be seen in Fig. 1. The main 
parameter that determines the number of hssion cycles is the rate of the r-process hux, 
which is a function of the location of the r-process path with respect to the stability 
line, and thus dictates the typical beta decay half lives. The location of the r-process 
path is a function of the neutron number density and temperature T, coming closer 
to stability for decreasing or increasing T, thus increasing cycle times and decreasing 
the number of cycles which in turn leads to smaller abundance predictions. 

Since the formation of heavy elements beyond Fe and Kr is a very sensitive measure of 
7], it can be used to provide an independent upper limit for the product Fig. 2 

shows observational (upper) limits to the primordial creation of heavy [7, 12, 13] and 
light [14, 15, 16, 17] element abundances to our results. The tightest constraints are 
given by the light elements including Li, Be, and B (see however recent doubts on the 
primordial ^Li abundance [18]) for which the conditions cannot differ much from the 
standard big bang. 

How do changes in the reaction rates leading to heavy elements affect our results? 
Some test calculations were performed with a variation of the ®Li(Q;,n)^^B rate, one of 
the two bottle necks toward heavy nuclei (the other one being ^^C(q;, 7)^®0). Recent 
experiments [19] seem to suggest that the rate used in our calculations [20] has to be 
increased by a factor of 3. It was found that such a change enters only linearly in 
the resulting heavy element abundances. This can be understood easily, as the seed 
production of heavy elements varies linearly with the Li-rate, and even for strong hssion 
cycling this behavior is not changed because Xheavy = ^seed x ■ A similar effect 
was found for the ^® 0 (n, 7)^®0 rate which was changed by a factor of 10 in a recent 
investigation [21]. Thus, the total change in heavy element abundances is a factor of 30. 
This is also shown in Fig. 2. (However, note that these changes are not included in the 
values given in Tab. 1). 

Provided that density huctuations exist with large scale lengths in comparison to the 
neutron diffusion length, the corresponding limits for 710 or llfehgQ change to 104 and 
1 . 6 , respectively, at which heavy element abundances are produced in inhomogeneous big 
bang models at a level comparable to the ones seen at lowest observable metallicities. 
This reduces the difference between the constraints from light and heavy elements, 
although the light element constraint is still tighter. It also underlines that not all 
reactions of importance are fully explored, yet, and future changes can be expected. 
Acknowledgement: TR is supported by the Alexander von Humboldt foundation. 
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Table 1: Mass fractions of heavy nuclei (see text). 
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Figure 1: r-process abundances for ryio = 416 and no diffusion (crosses) compared to the 
solar abundance (triangles). Fission cycling has enhanced the initial “seed” amount of heavy 
elements by 8 orders of magnitude. (The shift of the abundance peaks is due to the low neutron 
densities at late time [2,3]). 
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Figure 2: Limits on from light and heavy element abundances. Abundances are nor¬ 

malized to solar. Shown are the results for different length scales of inhomogeneities, i.e. case 
0 (full sq.), 1 (crosses), and 2 (open sq.). (The lines are merely drawn to guide the eye). 
Also shown is the result for the enhanced rates (dotted line). The horizontal solid line gives 
the observational upper limit for primordial heavy abundances. The limits on ^Ibh'^o resulting 
from the calculated values for the light element abundances are given by the vertical full and 
dashed lines. (See text). 
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